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Using Chiral Perturbation Theory, we investigate the volume dependence of the 
kaon semi-leptonic vector form factor at zero momentum transfer up to one loop 
order. We obtain the form factor for two volumes with L = 1.83 fm and L = 2.73 
fm and compare our results with the lattice data as quoted in Boyle et al. (Phys. 
Rev. Lett. 100 (2008) 141601). It turns out that our ChPT results in both volumes 
are in good accord with the lattice data such that the discrepancy is at most 1.1 
pecent. This study therefore, suggests that SU(3) ChPT is a reliable tool for lattice 
data, concerting vector form factor at zero momentum transfer, to extrapolate at 
large lattice size. 
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1 Introduction 



Lattice QCD calculations have improved significantly in recent years such that extraction of 
the parameters of the standard model are feasible, for one example see [H El El HI El E] . For 
a review on the lattice result on the low energy particle physics one may consult reference 
[7J. Among these parameters, the precision evaluation of the CKM matrix element, namely 
\V US \ is important from the vantage point of finding new physics footprint in the unitarity 
requirement of the first row of the CKM matrix [SI The main uncertainty in this 
unitarity relation is due to \V US \. From the measurements of the decay rate of the semi- 
leptonic kaon decay (K — > nev\ the so-called Ki 3 decay, one can only determine the 
combination \V us \f + (0), where the quantity /+(0) is the relevant vector form factor at zero 
momentum transfer. The average of measurements from different modes of Ki 3 is provided 
by PDG(2012)[I0] 

\V us \f + (0) = 0.21664 ± 0.00048 . (1) 

At our disposal are two approaches to determine \ V US \. On one hand, lattice QCD (LQCD) 
provides the vector form factor by numerical evaluation of the relevant functional integral 
of QCD, for instance in [4] it is found 

f+(0) LQC D = 0.9644 ± 0.0033 ± 0.0034 ± 0.0014, (2) 

where the first error is statistical and the second and third error are the size of the sys- 
tematic errors. As it is noted in [4], the uncertainty in the lattice data is dominated by 
the statistical, chiral and q 2 extrapolation. On the other hand, there is the application of 
an effective field theory, namely, chiral perturbation theory (ChPT). In this framework, 
the evaluation of f+(t) at one loop order is done by Gasser and Leutwyler [11]. However, 
the theoretical framework for the determination of the quantity V us is already discussed 
by Leutwyler and Roos in [12]. At one loop order, /+(0) is free from order p 4 low energy 
constants L\, hence, order p 4 corrections to the form factor purely arises from quantum 
corrections. There exist higher order ChPT based works [T3"l IT4"| ITS] [TB] which indicate the 
dominance of the two loop corrections to /+(0). It should be pointed out, however, that 
ChPT result at two loop order depends on the combination of the two order p 6 low energy 
constants C'[ 2 + Cg 4 . In conclusion, ChPT can provide model- independent prediction, but 
a precise value of /+(0) is not achievable due to uncertainty of the low energy constants. 

In this work we have made an attempt in order to make an estimate for the finite volume 
effects of /+(0) with non- vanishing spatial momentum in the framework of ChPT. These 
results can be considered as a guideline for lattice practitioners when they do the large 
volume extrapolation. In a previous work [17] , it was shown that the scalar form factor at 
the maximum momentum transfer acquires unexpected large finite volume correction for 
typical pion masses used in lattice calculations. The main motivation behind this work is 
actually to figure out the usefulness of ChPT application in finite volume for /+(0). In the 
literature can be found many works on employing ChPT to estimate the systematic errors. 
The finite size effects on the pion mass and pion decay constant in [IB"] [T9| [20] and on 
quark vacuum expectation values in [21] are examples for cases in which external momenta 
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are not involved. There is also a work for the meson matrix element in finite volume, with 
non-zero external spatial momentum in [22] . 

The rest of the article is organized as follows. In Sec. [2] chiral perturbation theory 
is briefly introduced. The next section introduces the strangeness-changing semi-leptonic 
kaon decay. The hadronic matrix element for the process in a tensor form at one loop 
order is recapitulated in Sec. HJ A short introduction to the application of ChPT in finite 
volume is given in Sec. |5]and all needed Feynman integrals in finite volume are calculated 
in Sec. [6j Finally, in Sec. Hour results for the vector form factor at zero momentum transfer 
are presented. We finish up with a conclusion. 



2 SU(3) chiral perturbation theory 

At low energies Quantum Chromo Dynamics (QCD) becomes a strongly coupled theory 
therefore the standard perturbation approach is no longer applicable. Chiral perturbation 
theory (ChPT) is an effective field theory to study the strong interactions at low energy. 
Spontaneously chiral symmetry breaking in QCD gives rise to Pseudo-Goldstone mesons 
which are considered as dynamical degrees of freedoms in the effective theory. ChPT is 
emerged in its modern in a paper by Weinberg [23J and developed latter on by Gasser and 
Leutwyler to higher order [2U |25] . External momentum, p 2 and quark masses, m q are the 
generic expansion parameters. At the lowest order SU(3) chiral Lagrangian contains two 
terms and takes on the form [23] 

£a = ^(«X + X+>, (3) 

where F Q is the pion decay constant at chiral limit and (...) = Trp(...) stands the trace 
over the flavors. We introduce the matrices w M and x± as following 

= iu^D^Uu^ — w , u 2 = U, 

x ± = u ] xu ] ± u\ ] u. (4) 



We can parameterize x i n terms of scalar and pseudo scalar external densities but in this 
work it is enough to set 



X = 2£? 



V 



\ 



m s J 



(5) 



The matrix U G SU(3) incorporates the octet of the light pseudo-scalar mesons 

U{(j>) = exp(zv / 20/F o ) , (6) 
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where 
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K° 



K+ \ 
K° 

2 m 



(7) 



We can obtain the relation m 2 = B (m u + from the lowest order Lagrangian. This 
relation allows us to take quark masses of order p 2 . In the covariant derivatives external 
fields are defined 



D u U = d u U-ir u U + iUL. 



(8) 



The left-handed and right-handed external fields are expressed by r M and Z M respectively. 
For the process we consider in this paper we just need to set 
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v: d w- 



(9) 



The weak coupling constant, g>2, is given in terms of Fermi constant and W mass by the 
relation g\ = Ay/^Gprn^. The next to leading order effective Lagrangian provided by 
[2U [25] contains twelve independent operators 

£ 4 = ii(«/) 2 + ^(n/KA,) + Ls^/m") + L 4 (wXXx+> 
+L 5 (uXx+> + L e(x+) 2 + MX-) 2 + ^(2L 8 + L 12 )(x+) 
+^(2L 8 - L 12 )( x 2 _) - iU{ft v u,u u ) + i(L 10 + 2L n )( WD 
-i(L 10 -2L 11 )</_^/r>, 



(10) 



where Li are the low energy constants which are obtainable phenomenologically and the 
field strength tensor is defined as 



F r M V ± v)F; 



R U, 



/r =u 

f£ v =dn v - d v F-i[F,r\, 



d"r v - <9V 



'IF 



3 The definition of the K — > n form factors 

Semileptonic weak decays of charge and neutral kaon known as are: 



(12) 
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K°(p)^n-(p')l + (p l )u l (p u ), (13) 

where subscribe / indicates electron or muon. There are two other processes which are the 
charge conjugate modes of the decays above. The matrix element for these processes, e.g. 
the neutral mode, 

K = |o ,, H(p f .p). ( 14 ) 

consists of two parts, keeping only the vector contributions. One part which defines the 
purely leptonic current 

r = u{p v )r{i-is)v{pi), (15) 

and the other part being our concern in this paper, incorporates the hadronic kaon-pion 
weak transition 

M,(p',p) = <7r-(p')\s^u(0)\K°( P )> . (16) 
The hadronic matrix element is generally defined by 



<n-(p')\siM0)\K°(p)> = ^[(p + pVfH0 + (p-P%/5 O7r ~(*)]- 



(17) 



A similar definition can be provided for the charge kaon. The two K i?) vector form factors 
f±° n (t) depend on the four momentum squared, t, transferred to the leptons 

t = {p-p') 2 = {pi+p v ) 2 . (18) 

The so-called scalar form factor as the S-wave projection of the matrix element can be 
defined as 

/oft) = /+(*)+ a 1 2 /-(*)• ( 19 ) 

m K ~ m i 

Given the definitions for the vector and scalar form factors, it is possible to obtain these 
dynamical low energy quantities in terms of temporal part and spatial part of the hadronic 
matrix element, namely, M-^ip' ,p)- For the vector form factor we find 

f , f) = (p i -p' i )Mo-(E p -E p ,)M t 

U V2(E piPi - ' 1 ) 

and for the scalar form factor we obtain 

f (f] _ f (m | I ( E m ~ EpPd (Pi + P'i)M - (E^+E f )Mj 

/0U /+U[ Mr - Ml (Epp'i - Ep Pi ) (pi - p'^Mo - (Ef - Ep)Mi 1 ' 

Where, Ai and Aii are respectively, the temporal and the spatial components of the weak 
vector current and Kaon and Pion energy are E$ = \Jm\ + p 1 and E^ = ^ M% + p* 2 , 
respectively. The relations above are useful when we look at the form factors in finite 
space in subsequent sections. 
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4 Analytical results for the weak matrix element 



The K13 form factors are known from [12J but in order to evaluate the form factors in finite 
volume one needs the matrix element in a tensor form. To this end in [T7j the matrix 
element at one loop order in the isospin limit is found 

M(p',p).e = ^[2g% + [^(m2) + ^(m2) + ^(m^)]r. e -[^(m^m^g 2 ^ 

+^ v (m 2 K , m% q 2 )]r u ^ + [-2{m 2 K - m\)U + ^{rn 2 K - m 2 )L 5 + -A(m 2 ) 
-J^ A ( m l) + ^Mjn 2 K )\q.e + B(ml,m 2 K , q 2 )(-^q 2 - j^m 2 K - ^ml)q.e 
+B(m 2 K ,m 2 v ,q 2 )(^q 2 - —m 2 K + ^m 2 )g.e - [^B^(m 2 n ,m 2 K ,q 2 ) 
+ ^£vKo m 2 , q 2 )}q»e u + ^(m 2 , m 2 K , q 2 ) [^(p + p'Yq.t + ^Q-e 
+ l m l^ + \ m l^ ~ g^l + B ^m 2 K , ml, q 2 )[^(p + p'fq.e 
+\q»q.e + -m 2 K e» - -m 2 ^ - -gV]l , (22) 



where 



6 A 6 n 2 



r = p' + p , q = p — p' , (23) 



and e is the polarization four vector of the W boson. In the expression above scalar 
integrals A and B as well as tensor integrals B^ and B^ v are introduced in the Appendix. 
The temporal and spacial parts of the matrix element, A4^(p',p), are then available from 
Eq. 1221 These are needed in our evaluation of the form factors in finite volume as defined 
in Eq. [20] and Eq. [2TJ The relevant Feynman integrals in finite volume are calculated in 
section |6j 

5 ChPT application in finite lattice box 

In lattice QCD, simulations are performed in a cubic volume (V=L 3 ) with periodic bound- 
ary conditions imposed on the hadronic fields 

(f)(x) = (f)(x + nL) . (24) 

Therefore the three vector momenta of hadrons become discrete 

p=—n, (25) 

where, n is a three dimensional vector with integer components (n x , n y , n z ). The application 
of chiral perturbation theory to study finite volume effects are introduced in original works 
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by Gasser and Leutwyler, see [26J, [27J I2H] for detailed discussion in this regard. It is 
important to note that with the periodic boundary condition, the effective Lagrangian in 
finite volume is the same as the one in infinite volume. Finite volume corrections get their 
effects from modification of the hadron propagation in space-time. Given the quantization 
of the momenta in finite volume, the two-point correlation function becomes 



Our power counting quantity in finite volume calculations is the quantity m^L where it 
turns out that the zero mode of the pion filed is not strongly coupled if the condition 
m n L » 1 is fulfilled. This is the so-called p-regime. In addition, ChPT gives reliable 
results when F n L » 1. 

6 Feynman integrals in finite volume 

In this part we calculate the Feynman integrals in finite volume for a generic momentum 
transfer q 2 . There are two types of integrals which appear in our expressions for form 
factors: scalar integrals and tensor integrals. In fact we wish to evaluate finite volume 
correction for a given integral for which we define AI = Iy — 1^, where subscripts oo and 
V indicate integration in infinite and finite volume respectively. 

6.1 One loop scalar integrals 

The simplest integral we encounter in this work is related to the tadpole Feynman diagram 



In finite volume, momentum gets quantized and therefore integration over momentum is 
replaced with summation over momentum 




(26) 



p 




(27) 





In obtaining the second line, we have employed the Poisson summation formula 
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(29) 



6 



11 


1 


2 


3 


4 


5 


6 


7 8 


9 


10 


11 


12 


13 


14 


15 


16 


17 


18 


19 


20 


m(n) 


6 


12 


8 


6 


24 


24 


12 


30 


24 


24 


8 


24 


48 





6 


48 


36 


24 


24 



Table 1: The multiplicity factors m(n) are provided for 1 < n < 20. 



By taking a contour integration over p and then performing the three dimensional integral 
we achieve the known result 



AA = ~4^L £ T^-M^MMLln]), (30) 

n=l ' ' 

where, K\ is the modified Bessel function of order one and the multiplicity factor m(n) 
stands for the number of possibilities that the relation n = n\ + n 2 + n| is satisfied for a 
given value of n with positive and negative integer numbers of ni, n 2 and n^. m(n) factors 
are listed in Table. [TJ 

The next integral we should evaluate in finite volume as a new one is related to the 
rescattering effects at momentum transfer q 2 

B(m 2 ,M 2 ,q 2 ) = - [ 7^7-^ 7-77^ ^ — (31) 

v ' lH ' i J (27v) d (p 2 -m 2 )((q + p) 2 - M 2 )' v ; 

where we assume hereafter that M > m. With the application of the Poisson summation 
formula we have 

i ^ f dp 1 



B v (m*,M 2 ,q 2 ) = -^j:J 

v 



2tt (p 2 -m 2 )((p + q) 2 - M 2 ) 



r dp ^ f d 3 p i e iLp.n 



2vr ^ J (2n)3(p 2 -m 2 )((p + q) 2 -M 2 y 

(32) 

Making use of the Feynman parameter formula in followed by redefining the variable po 
we will arrive at 



A£?(™ 2 , M 2 , q 2 ) = -i /* dx [ ^ Y 

Jo J 2tt^ 



d 3 p 



e 



iLp.n 



(27r) 3 [pi — (p+ (1 — x)q ) 2 + x(l — x)q 2 — xm 2 — (1 — x)M 2 ] 2 



.(33) 



At the next step, we begin by taking the contour integral over p and then make a redefi- 
nition of the variable p to obtain 

1 rl r ^j? iLp.n 

AB(m 2 ,M 2 ,q 2 ) = - W dx e"^" / £L— ^— rvr^m 

4 H^o {27r) d \p z — x{l — x)q z + xm z + (1 — xjM^l 6 ' 2 
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11 


C n (ax) 


n 


C n (ax) 


1 


2 Cos(ax) + 4 


11 


8 Cos(3ax) + 16 cos(ax) 


2 


8 Cos(ax) + 4 


12 


8 Cos(2ax) 


3 


8 Cos{ax) 


13 


8 Cos{3ax) + 8 Cos(2ax) + 8 


4 


2 Cos(2ax) +4 


14 


16 Cos{3ax) + 16 Cos{2ax) + 16 Cos(ax) 


5 


8 Cos (2ax) + 8 Cos(ax) + 8 


15 





6 


8 Cos(2ax) + 16 Cos(ax) 


16 


2 Cos(4aa;) + 4 


7 





17 


8 Cos(4ax) + 8 Cos(3aa;)+ 
16 Cos(2ax) + 8 Cos(aa;) + 8 


8 


8 Cos(2ax) + 4 


18 


6 Cos(4ax) + 10 Cos{iax) 
+ 16 Cos(ax) +4 


9 


2 Cos(3ax) + 16 Cos(2ax) + 8 Cos(ax) + 4 


19 


18 Cos(3ax) + 6 Cos(ax) 


10 


8 Cos(3ax) + 8 Cos (ax) + 8 


20 


8 Cos(Aax) + 8 Cos(2ax) + 8 



Table 2: Functions C n (aa;) are provided for 1 < n < 20 when external momentum g = 
g x (1,0,0) is chosen. 



The exponential factor e~ tLx ^ n explicitly breaks the rotational symmetry in the expression 
above. We carry out the integral over the vector momentum in two steps. We take first 
an integral over the angular part of the three dimensional momentum and then we make 
use of the convolution technique to perform the final integral. We find the following result 

A£?(™ 2 , M 2 , q 2 ) = -L f 1 dx C n (ax) K (w Q) , (35) 

where w = L\n\ and Q = \Jxm 2 + (1 — x)M 2 — x(l — x)q 2 . K is the modified Bessel 
function of rank one. This is a generalization of the case with M = m obtained in |29j . 
Functions C n {ax) introduced in the expression above involve the exponential factor e - lLx i- n 
where we have summed over all possible ways that for a given n the relation n = nf + n^+n 2 
is satisfied. In this article we take the the momentum transfer, q, along the x- axis, i.e, 
q = (^,0,0). We provide functions C n (ax) in Tabled where a = Lq x . It is easy to see 
that for a = 0, functions C n (ax) are identical to multiplicity factors m{n). 



6.2 One loop tensor integrals 

In this part we consider the tensor integrals by calculating the temporal and spatial com- 
ponents. We begin with the temporal component of the tensor integrals. The next integral 
we then need to consider in finite volume is 

B° is the temporal component of the tensor integral defined in Appendix A. By repeat- 
ing the procedures stated above we can readily prove that 

AB°(m 2 , M 2 , ? 2 ) = -^Ef&x C n (ax) K (w Q) , (37) 
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11 


D n (ax) 


n 


D n (ax) 


1 


2 Sin(pt.x) 


11 


8 5m(3ax) + 16 Sin(ax) 


2 


8 Sin(ax) 


12 


8 »Sm(2ax) 


3 


8 Sin(ax) 


13 


8 Sin(3ax) + 8 Sm(2aa;) 


4 


2 Sin(2ax) 


14 


16 Sin(3ax) + 16 Sin(2ax) + 16 Sin(ax) 


5 


8 S f m(2o;a;) + 8 Sin(ax) 


15 





6 


8 i5in(2aar) + 16 Sm(aa;) 


16 


2 S , m(4ci!a;) 


7 





17 


8 5m(4aa;) + 8 Sin(3ax) 
+16 S , m(2aa;) + 8 Sin(ax) 


8 


8 Sin(2ax) 


18 


6 S'm(4Q!a;) + 10 Sin(3ax) + 16 5m(aa;) 


9 


2 S'm(3aa;) + 16 Sin(2ax) + 6 Sin(ax) 


19 


18 Sin(3ax) + 6 Sm(air) 


10 


8 5m(3aa;) + 8 Sin(ax) 


20 


8 Sin(Aax) + 8 S'zn(2aa:) 



Table 3: Functions D n (ax) are provided for 1 < n < 20 when external momentum g* = 
q x (1,0,0) is chosen. 



where, go — E p i — E p . Now we look at the temporal component of the tensor integral 

B«V, M 2 , g 2 ) = 1 / ( 2^ (p 2_ m2)( J° +p)2 _ M2) - ( 38 ) 
We repeat the procedures sketched above and finally arrive at 

AB 00 (m 2 , M 2 , g 2 ) = -— L Enf^ Q ^ 

8vr 2 L ^ |n| Jo 

^ |4| f dx x 2 C n (ax) K (w Q). (39) 



Moreover, the tensor integral has a spatial component 

*(Mt) = \!^ V-rn*)£+pr-WV (40) 

In order to find the integral B x in finite volume we follow the same path as we did to 
evaluate B°. Our final result read 

A^(m 2 ,M 2 ,g 2 ) = -L £ f 1 dx Q K x {w Q) 

o7r 2 r^ Q Jo \n\ 

[ x dx x C n (ax) K (w Q). (41) 



Qx 



n^0 ' 



Functions D n (ax) appear when for a given n, we compute the summation J2i n x e~ lLxq ' n 
over all possible values of (n^na, TI3) that fulfill the relation n = n 2 + n 2 , + n\. Functions 
D n (a) are listed in Table El The last integral in finite volume is the spatial component 
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11 


F n (ax) 


n 


F n (ax) 


1 


2 Cos(ax) 


11 


72 Cos(iax) + 16 Cos(ax) 


2 


8 Cos(ax) 


12 


32 Cos{2ax) 


3 


8 Cos(ax) 


13 


72 Cos(3a:c) + 32 Cos(2ax) 


4 


8 Cos(2ax) 


14 


144 Cos(3ax) + 64 Cos(2ax) + 16 Cos(aa;) 


5 


32 Cos(2aa;) + 8 Cos(ax) 


15 





6 


32 Cos(2ax) + 16 Cos(ax) 


16 


32 Cos{4ax) 


7 





17 


128 C*os(4ax) + 72 Cos(3ax) 
+64 Cos (2ax) + 8 Cos(ax) 


8 


32 C*os(2aa;) 


18 


96 Cos(4aa;) + 90 Cos(3ax) + 16 Cos(ax) 


9 


18 Cos(3ax) + 64 Cos(2ax) + 6 Cos{ax) 


19 


162 Cos(3aa;) + 6 Cos(ax) 


10 


72 C*os(3a:z;) + 8 Cos(ax) 


20 


128 C*os(4ax) + 32 Cos(2ax) 



Table 4: Functions -F n (ax) are provided for 1 < n < 20 when external momentum q = 
q x (1,0,0) is chosen. 



B xx of the tensor integral defined in the Appendix 

5-(m 2 , M 2 , g 2 ) = - / „ -r . (42) 



z7 (27r) d (p 2 -m 2 )((g + p) 2 -M 2 )' 
For the finite volume correction of the integral above we obtain 

A5-(m 2 ,M 2 ,g 2 ) = -L £ f dx [-^M Q Q) 

Fn(ax) -Q 2 [*o(«> Q) + # 2 (™ Q)] + Mw q)\ 



2\n\ 2 " " L\n\ 



ql 
4 ' 



8tt 2 - „ 

q x 



An 2 , n 



Y") / x 2 <ix C„(ax) -K"o(w Q) 

Y t xdx Q Kt(w Q), (43) 



Functions F n (ax) listed in Table. HI are obtained by computing J2 n x e tLxq ' n for a given n 
over all possible ways that the relation n = nf + nf + n| holds. 



7 /+(0) in finite volume 

At this point we present our numerical results. In the calculations involving /+(0) in 
infinite volume, we use as input parameter for the pion decay constant 

F T = 0.0924 GeV. (44) 

Since we look at the vector form factor at zero momentum transfer, the low energy constant 
Lg does not show up in our expression as it is evident from Eq. [221 We show in Fig. [1] the 



10 



itl = 200 MeV, m K = 500 MeV — 

m„ = 300 MeV, m K = 500 MeV 

m. = 400 MeV, m K = 500 MeV 



Figure 1: The ratio Rf + (o) is plotted versus the spatial size of the volume, L, for three 
different values of pion masses with fixed kaon mass. 



results for the ratio Rf + (o), defined as 



R 



/+(o) 



■ff(0)-/j°(0) 



(45) 



in terms of the linear size of the volume, L. It is worth to note that the ratio has no 
dependence on Lg. In the figure different lines stand for different choices for the pion 
mass i.e m n = 0.2 GeV, = 0.3 GeV and m n = 0.4 GeV with the fixed value for kaon 
mass, ttik = 0.5 GeV, while for the eta mass we use the GMO relation at leading order, 
m 2 = (Am 2 K — m^)/3. They indicate two standard characteristic features since the ratio 
tends to zero asymptotically and on top of that the ratio grows with decreasing the pion 
mass. 

Moreover, we have calculated the vector form factor at zero momentum transfer for 
two ensembles corresponding to volumes with L = 1.83 fm and L = 2.74 fm as quoted in 
[4J. For the smaller volume with L = 1.83 fm, in Table. Infinite volume corrections of the 
vector form factor are presented for both cases when eta-kaon loops are ignored and when 
both pion- kaon and eta-kaon loops are present. Our finite volume corrections are in good 
accord with the lattice results when eta-kaon loops are subtracted and agree a little less 
when eta-kaon loops are included. The size of the discrepancy changes from 0.8 percent 
for the lightest pion mass i.e. m n = 0.428 GeV to 0.25 percent for the heaviest pion mass 
i.e. m n = 0.674 GeV. It is an important observation that lattice data are accurate enough 
to distinguish the effects of the eta propagation in finite volume. 

On the other hand, the results in Table. El corresponding to volume with L = 2.74 fm, 
also show a good agreement with the lattice data, specially for the heavier pion masses. 
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JT(o) 


/|(0) 


/+(0) 


A/+(0) 


A/+(0) 


mo) 


(GeV) 


(Pole) 0] 


(Quadratic) [I] 


(ChPT) 


Without r\ loops 


With r\ loops 


With r\ loops 


0.428 


0.9889 (26) 


0.9866 (33) 


0.99315 


-0.01270 


-0.01447 


0.97868 


0.557 


0.9951 (6) 


0.9959 (9) 


0.99760 


-0.00497 


-0.00608 


0.99152 


0.674 


0.99925 (8) 


0.99938 (12) 


0.99950 


-0.00215 


-0.00278 


0.99672 



Table 5: The vector form factor for different pion masses in a finite volume with spatial 
size L = 1.83 [fm] in the last column is compared with the same quantity in the third and 
forth columns, taken from [I], evaluated within lattice QCD using pole and quadratic fit, 
respectively. 



m n 
(GeV) 


mo) 

(Pole) g] 


(Quadratic) [4] 


/+(0) 
(ChPT) 


A/+(0) 
without r/ loops 


A/+(0) 
with r\ loops 


_ /|(0) 
with r) loops 


0.329 


0.9774 (35) 


0.9749 (59) 


0.98825 


-0.00233 


-0.00241 


0.98583 


0.416 


0.9841 (29) 


0.9806 (39) 


0.99301 


-0.00086 


-0.00093 


0.99207 


0.556 


0.9960 (7) 


0.9962 (9) 


0.99769 


-0.00018 


-0.00022 


0.99747 


0.671 


0.9991 (2) 


0.9990 (2) 


0.99952 


-0.00005 


-0.00006 


0.99946 



Table 6: The vector form factor for different pion masses in a finite volume with spatial 
size L = 2.74 [fm] in the last column is compared with the same quantity in the third and 
forth columns, taken from [I], evaluated within lattice QCD using pole and quadratic fit, 
respectively 



The result indicate a deviation of about 1.1 percent for lightest pion mass, i.e. = 0.329 
GeV down to 0.04 percent for the heaviest pion mass, i.e. m n = 0.671 GeV. The eta loop 
contribution is again very small, but its effect can be traced in the lattice data. 

8 Conclusion 

In this work we have found the kaon semi-leptonic vector form factor in a finite box at 
zero momentum transfer. Comparing with the available lattice result, we understand that 
the application of the SU(3) ChPT in finite volume for the discussed process is successful. 
In fact, in both volumes for the larger pion mass our results and lattice data differ by less 
than 0.5 percent. In addition, we realize that the result for the scalar form factor in finite 
volume at zero momentum transfer obtained in this article is much more satisfactory than 
the scalar form factor at the maximum momentum transfer as found in |17j . 

9 Appendix 

We introduce the necessary one loop scalar Feynman integrals for the discussed decay 

. , is 1 f d d P 1 , x 

Aim 2 ) = - / t-^-j-, r, 46 

v ' i J (2vr)<V -m 2 ' v ' 

12 



B(m 2 ,M 2 ,q 2 ) =-/t^ ^ — ^- ^ ^ (47) 



(27r) d (p 2 - 


m 2 )((p + g) 2 — 


M 2 )' 


d d p 






(27l) d (p 2 - 


• m 2 )((p + g) 2 - 


•M 2 )' 









and for the tensor Feynman integrals 



B^(m 2 ,M 2 ,q 2 ) -J (2 ^ d (p2 _ m 2)(^7 g) 2 - M 2 ) ' ( 49 ) 
One can write the tensor integrals in terms of scalar functions by applying lorentz symmetry 

B^m 2 ,M 2 ,q 2 ) = q.B^m 2 , M 2 ,q 2 ), (50) 

S^(m 2 , M 2 , q 2 ) = q^q u B 12 (m 2 } M 2 , g 2 ) + g^B 22 (m 2 , M 2 , g 2 ). (51) 
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